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We report a detailed investigation of RECoPO (RE = La, Pr) and LaCoAsO materials performed 
by means of muon spin spectroscopy. Zero-field measurements show that the electrons localized on 
the Pr 3+ ions do not play any role in the static magnetic properties of the compounds. Magnetism 
at the local level is indeed fully dominated by the weakly-itinerant ferromagnetism from the Co 
sublattice only. The increase of the chemical pressure triggered by the different ionic radii of La 3+ 
and Pr 3+ , on the other hand, plays a crucial role in enhancing the value of the magnetic critical 
temperature and can be mimicked by the application of external hydrostatic pressure up to 24 kbar. 
A sharp discontinuity in the local magnetic field at the muon site in LaCoPO at around 5 kbar 
suggests a sizeable modification in the band structure of the material upon increasing pressure. 
This scenario is qualitatively supported by ab-initio density-functional theory calculations. 

PACS numbers: 74.70.Xa, 71.15.Ap, 75.50.Cc, 76.75. +i 



I. INTRODUCTION 

The discovery of high-T c superconductivity in F-doped 
REFeAsO (1111) was one of the most significant break- 
throughs in the fiel d o f condensed matter physics dur- 
ing the last decades A striking worldwide interest has 
been indeed recently devoted to the synthesis of this class 
of materials, looking for suitable chemical substitutions 
of rare-earth (RE) ions, transition metals (TM) and pnic- 
togen elements (Pn) allowing to enhance the critical tem- 
perature T c , the highest value T c ~ 55 K being currently 
reached in SmFeAsOi_ x F x pfil Electron doping can be 
realized directly on the FeAs layers also leading in turn 
to superconductivity. In this respect, one of the most 
studied chemical substitution of TM is Fei_ K Co x . This 
applies to Fe-based materials belonging both to 1111 and 
to 122 families, the parent compound for the latter case 
being BaFe 2 As 2 P tn l 

Beyond superconductivity, 1111 materials show in- 
teresting magnetic features associated with the mutual 
interaction among localized electrons onto the exter- 
nal shells of RE ions and itinerant carriers from the 
TM sublattice.^ A strong f — d hybridization was 
shown to be present, for instance, in superconducting 
SmFeAsOi-^Fa; under conditions of optimal doping! 1 - 4 ^ 
In the case of undoped REFeAsO compounds, a spin 
density wave (SDW) phase below T N ~ 140 K associ- 
ated with itinerant d electrons from Fe is found to coex- 
ist with an antiferromagnetically (AFM) ordered phase 
of RE magnetic mome nts at much lower temperatures 
(TRE ~ 5 - 10 K)PE1 it should be stressed how the full 



substitution of Fe by Co has attracted particular interest 
e. g. in RECoAsO compounds where much more com- 
plex magnetic behaviours were shown by means of both 
macroscopic and local experimental techniques EHUD xhe 
Co sublattice, in particular, is known to enter a weakly- 
itinerant ferromagnetic (FM) phase below Tq — 60 — 80 
K, the precise v alue of Tc being strongly dependent on 
RE) 17 * 19 ! 20 * 25 ! 26 ! At lower temperatures, the occurrence 
of other phase transitions is clearly observed in sam- 
ples containing magnetic RE ions as a consequence of 
the strong interaction among the two magnetic sublat- 
tices. This behaviour, in fact, is typically interpreted as 
the result of progressive FM-AFM transitions of the Co 
sublattice induced by the RE ions followed at the low- 
est tempe ratures by the full magnetic ordering of the RE 
sublatticePEEHSI 

A similar phenomenology has been recently r eported 
in the isostructural P-based RECoPO compounds! 17 * 27 * ^ 
LaCoPO has been investigated by means of 31 P nuclear 
magnetic resonance (NMR) showing that Moriya's the- 
ory of self-consistently renormalized spin fluctuations for 
weakly-itinerant magnets (Tc ~ 35 K) well describes the 
experimental results!22H2U Magneto-transport, dc magne- 
tometry and NMR measurements on compounds with 
magnetic RE ions like Sm 3+ and Nd 3+ show the occur- 
rence of multiple mag netic p hase transitions similarly to 
the case of RECoAsO PEl i n RECoPO compounds an 
interaction between the two magnetic sublattices much 
stronger than what is occurring in RECoAsO should be 
expected. The P/As isovalent substitution for the Pn el- 
ement, in fact, is known to introduce a strong chemical 
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pressure making the RE ions much closer to the itiner- 
ant layers. In the case of CeFeAsi-^P^O, for instance, 
these effects are known to gradually suppress the SDW 
phase associated with Fe upon increasing x, driving at 
the same time the RE sublattice from an AFM ordered 
phase through a FM ground state (GS) and finally to- 
wards a Kondo-screened phase where heavy- fermion phe- 
nomenology is recovered pHEU 

In order to evidence how crucial the role of the chemical 
pressure on the magnetic properties of RECoPO mate- 
rials is, we performed measurements of muon spin spec- 
troscopy on these compounds. In this paper we report on 
the results obtained in LaCoPO and PrCoPO, where the 
local magnetism was investigated upon the application 
of external hydrostatic pressure up to 24 kbar. Remark- 
ably, in spite of the high value of the magnetic moment 
expected for the free Pr 3+ ion (/ip r ~ 3.6 /ub), electrons 
localized on the external shells of Pr 3+ do not play any 
role in the local static magnetic properties of PrCoPO 
probed by muons. For both compounds, only the itin- 
erant ferromagnetism from the Co bands dominates the 
observed response. The chemical pressure triggered by 
the full Pr 3+ /La 3+ substitution, on the other hand, has 
a sizeable effect in enhancing the value of the critical 
temperature of the itinerant ferromagnetic phase. In 
this respect, external hydrostatic pressure is shown to 
lead to the same result as Pr/La substitution both in 
LaCoPO and, to a lesser extent, in the isostructural com- 
pound LaCoAsO. Furthermore, both chemical and exter- 
nal pressures strongly suppress the local magnetic field 
at the muon site while leaving the magnetic moment per 
Co ion substantially unchanged. It should be remarked 
that the chemical shrinkage of the lattice is intrinsically 
expected to be characterized by a higher degree of non- 
hydrostaticity that, in the case of 122 systems, was shown 
to play a drastic role in governing the resulting magnetic 
properties.^ However, claims of close analogies among 
the effect of chemical and external pressures were re- 
ported concerning both 1111 and 122 compoundsPsHSni 
The results presented in this paper further confirm this 
latter scenario also for Co-based 1111 materials. Ab ini- 
tio density-functional theory calculations have been de- 
veloped to describe the effects of pressure both on the in- 
terstitial crystallographic sites for the muons and on the 
electronic bands of LaCoPO. Results support the exper- 
imental findings and suggests that chemical and external 
pressures can both trigger a change in the electronic band 
structure. As a result, chemical and external hydrostatic 
pressures act similarly on the magnetic properties of La- 
CoPO. 



II. EXPERIMENTAL DETAILS 

Loose powders of LaCoAsO, LaCoPO and PrCoPO 
were grown via solid-state reactions as described in detail 
in Refs. l2~4l and |2"51 The structural properties of the lat- 
tices were measured at room temperature by means of a 
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FIG. 1: (Color online) Observed (red circles) and calculated 
(blue solid lines) X-ray powder-diffraction patterns at room 
temperature for the investigated samples of LaCoAsO (see (a) 
panel), LaCoPO (see (b) panel) and PrCoPO (see (c) panel). 
Black lines are best-fits to experimental data according to a 
Rietveld analysis. 



Rigaku X-ray diffractometer with Cu K Q radiation. The 
Rietveld analysis of the diffraction patterns (see Fig. [IJ 
confirmed that all the samples crystallized in the tetrag- 
onal phase, space group P4/nmm, and allowed us to ex- 
tract the values of the lattice parameters a and c (see 
Tab. [[]). It is clear from Tab. |l]how both the P/As and 
Pr/La substitutions increase the chemical shrinkage of 
the cell and reduce both a and c accordingly. This is 
in agreement with what is genera lly reported for all the 
1111 family of compounds! 19 * 34 * ^ 

Measurements of zero-magnetic-field (ZF) muon spin 
spectroscopy (/x + SR) were performed at the GPD spec- 
trometer (/iEl beamline) of the S/iS muon source at the 
Paul Scherrer Institut, Switzerland (see Ref. H21 for a 
comprehensive introduction to /i + SR techniques). Pres- 
sure (P) was applied at ambient temperature (T) in a 
double-wall piston-cylinder cell made of MP35N alloy 
and its value was quantified by ac susceptometry at T ~ 3 
K from the shift of the superconducting critical temper- 
ature of a small In wire inside the cell. The transmitting 
medium Daphne oil 7373 was employed assuring that the 



TABLE I: Lattice parameters for the investigated samples af- 
ter Rietveld refinements of X-ray powder diffraction patterns 
displayed in Fig. I] 



Compound 


a (A) 


C (A) 


Volume (A 3 ) 


x 2 


LaCoAsO 


4.048(8) 


8.462(7) 


138.73(1) 


3.32 


LaCoPO 


3.966(7) 


8.365(0) 


131.62(7) 


2.12 


PrCoPO 


3.9208(4) 


8.212(4) 


126.25(0) 


2.31 
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pressure conditions were always nearly hydrostatic in the 
experimental rangeP^H ^he maximum P value attain- 
able at low T with this setup is close to 24 kbar. The 
pressure cell (PC) intrinsically leads to an extremely high 
background level in ZF-/i + SR measurements, whose be- 
haviour as a function of T was tested and measured in 
an independent set of experiments. For this reason, a 
ZF-/i + SR characterization of the samples was also per- 
formed at the low-background spectrometers Dolly and 
GPS (7rEl and 7rM3 beamlines, respectively) at PSI. 

dc magnetization measurements were performed by us- 
ing the commercial superconducting quantum interfer- 
ence device (SQUID) magnetometer MPMS-XL7 (Quan- 
tum Design). A piston-cylinder CuBe PC (EasyLab 
Mcell 10) was used to apply hydrostatic P < 11 kbar. 
Again, the Daphne oil 7373 was employed as transmit- 
ting medium. P was applied at ambient T and its value 
quantified at low T from the shift of the superconducting 
critical temperature of a small Sn wire inside the cell. 



III. EXPERIMENTAL RESULTS 

The measured spin depolarization functions for the im- 
planted muons (/i + ) as a function of time (t) in ZF con- 
ditions were fitted for all the samples and at all the in- 
vestigated T values by the general expression 



A T (t) = A 



+ (l-a PC )G s T (t) 



(1) 



Here the amplitude etpc accounts for the fraction of 
incoming /i + stopping in the PC. The depolarization 
of these /i + is a Gaussian-like relaxation governed by 
a nearly T-independent ape term arising from nuclear 
magnetism inside the MP35N alloy. The remaining frac- 
tion (1 — ape) of fj, + is implanted directly into the sample 
and, accordingly, the relative depolarization G|*(i) can be 
described as 



G s T (t) = [l-U m (T)]e" 



+ 



a 1 - (T) J F(t)£i- L (t) + J (T)T>" (t) 



(2) 



Here the quantity V m (T) represents the fraction of fi + ex- 
periencing a static local magnetic field or, equivalently, 
the magnetic volume fraction of the investigated sam- 
ple. In the paramagnetic limit, namely V m (T) — 0, no 
static field of electronic origin contributes to the depolar- 
ization and only the weak contribution from the nuclear 
magnetic moments leads to a slow Gaussian depolariza- 
tion with characteristic rate a>j (typical measured values 
(Tn ~ 0.1 /is -1 ). Below the magnetic-order critical tran- 
sition temperature Tq, the superscript _L (||) refers to 
fi + experiencing a local static magnetic field in a per- 
pendicular (parallel) direction with respect to the initial 
ji + spin polarization. The amplitudes a ±: ^ must then 
satisfy the condition [a 1 - (T) + a" (T)] = V m (T) accord- 
ingly. In the presence of a long-range magnetic order 



inside the sample, a coherent precession of /i + around 
the local magnetic field can be discerned in the a 1 - 
amplitude and described by the oscillating function F(t). 
The (either exponential or Gaussian) damping function 
D 1 - (t) reflects a distribution of local magnetic field values 
at the /i + site. The a" component, on the other hand, 
is typically damped by the exponentially-decaying func- 
tion £>"(£) = e~ A "* probing spin-lattice-like relaxation 
processes (A" ~ 0.1 ^ s *)• 

Standard oscillating functions F(t) = cos (jB^t + (f>), 
where 7 = 2ir x 135.54 MHz/T is the magnetogyric ratio 
for are in good agreement with experimental data 
(statistical x 2 ~ 1 — 1.2). Negative values for the phase 
1 0| ~ 20° — 30° were systematically measured at the 
low-background spectrometer Dolly, similarly to what re- 
ported in literature for RECoAsO compounds. 26 Phase 
values close to —30° are generally accepted to be an ex- 
perimental evidence of magnetically-ordered phases in- 
commensurate with the underlying crystalline lattice. 

The three main contributions to the local field at the 
muon site B^ in the case of ferromagnetic materials come 
from the dipolar field, the transferred hyperfine field and 
the Lorentz field As it will be discussed in more de- 
tail later in Sect. IV (see Eq. Q), the crucial physical 
parameters governing the amount of those contributions 
are the distance between the /x + and the magnetic ions, 
the magnetic moment of the ordered magnetic phase and 
the density of spins at the fi + site.^ Once the interstitial 
crystallographic position of the \x + is known, the inde- 
pendent knowledge of the value of the magnetic moment 
(by means of e. g. dc magnetometry) is of crucial impor- 
tance since it allows one to directly access the transferred 
hyperfine fieldpSSO 

Representative ZF /i + -depolarizations for LaCoPO ob- 
tained at Dolly (ambient P) are presented in the inset of 
Fig. [2j Coherent oscillations indicative of the estab- 
lishment of a long-range magnetic order can be clearly 
discerned for T < Tq — 33 K down to the lowest in- 
vestigated temperature. Experimental data were fitted 
by means of Eq. ([I]) accordingly, where ape = 0. The 
T-dependence of the magnetic volume fraction V m (not 
shown) confirms that the sample is fully magnetic below 
To The B^ vs. T trend is plotted in the main panel 
of Fig. [2] and can be fitted according to a power-law 
function 



B^T) = B„(0) 1 - 



T 

Tc, 



(3) 



eventually saturating at B^(Q) = 587.5 ± 2.5 G, the ex- 
ponent turning out to be j3 = 0.34 ± 0.01. 

Also in the case of PrCoPO, measurements performed 
at Dolly reveal a static long-range magnetism probed by 
/i + in 100% of the magnetic volume (not shown) below 
Tq ~ 48 K. The value for Tc is sizeably increased with 
respect to what is observed in LaCoPO, consistently with 
what is reported for RECoAsO compoundsP3 Remark- 
ably, the B^ vs. T trend still can be well-described by 
Eq. ([3]) where again ft — 0.34 ± 0.01 (see the main panel 
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FIG. 2: (Color online) Main experimental results of ZF-/x + SR 
(ambient P, Dolly spectrometer). Main panel: £? M vs. T for 
LaCoPO and PrCoPO. Dashed lines are best-fits to experi- 
mental data according to Eq. [3] Inset: raw ZF spin depolar- 
ization data for LaCoPO at some selected T values. Contin- 
uous lines are best fits to experimental data according to Eq. 
[l] where ape = 0. 



FIG. 3: (Color online) Main panel: B M (T) at different P 
for LaCoPO. Dashed lines are best fits to experimental data 
according to Eq. [3] with (3 = 0.34 as fixed parameter. Close 
and open symbols for ambient P data refer to measurements 
performed at Dolly and at GPD with unloaded cell, respec- 
tively. Inset: Tc vs. P. The continuous line is a best-fit to 
data according to a linear function. 



of Fig. |2j. In the case of PrCoPO, only the local field 
Bfj,{0) = 132 ± 2 G is more than 4 times lower than in 
LaCoPO. In spite of the high value for the localized mag- 
netic moment expected for the free Pr 3+ ion (~ 3.6 
the results for the PrCoPO sample are qualitatively iden- 
tical to the case of LaCoPO. This is clearly different from 
what was previously reported in the literature for the 
As-based compound PrCoAsO, where a more complex 
phenomenology for ZF-/Lt + SR (displaying a splitting of 
internal fields) was observed.^ Nevertheless, dc magne- 
tization measurements for LaCoAsO and PrCoAsO were 
reported to display no sizeable differences amon g th em- 
selves, no anomalies being detected for T < To The 
present results (together with what is reported in Ref. 
|2"B")) then suggest that both in PrCoPO and in PrCoAsO 
the magnetic moments localized on the Pr 3+ ions do not 
play any role in the static magnetic features of these com- 
pounds. The phenomenology observed in PrCoAsCPsl 
then can possibly be due to the presence of two different 
/x + sites whose statistical population is strongly depen- 
dent on T in the investigated experimental range. 

The overall scenario deduced from our measurements 
on LaCoPO and PrCoPO and from the previous char- 
acterization of the materials (see Sect. |n| suggests that 
the main effect of the full Pr/La substitution is only as- 
sociated with a lattice shrinkage induced by the different 
chemical pressures associated to the different ionic radii 
of La 3+ and Pr 3+ . This is different from the cases of 
all other magnetic RE ions where the interplay of the 
localized and the itinerant magnetic degrees of freedom 
leads to progressive FM-AFM re-ordering effects on the 



Co sublattice. Therefore, it seems safe to deduce that 
p r 3+ ma g ne tic moments do not contribute to the local 
static magnetic properties of the compound that are only 
governed by weakly-itinerant FM from Co. 

The validity of these arguments is corroborated by the 
measurements under external P performed at the GPD 
spectrometer. In spite of the fact that almost half of the 
incoming /i + stop inside the pressure cell (ape ~ 0.5), 
clear oscillations can be distinguished from the ZF sig- 
nal associated with all the samples below characteristic 
critical transition temperatures Tq(P). The analysis of 
V m (not shown) clearly shows that the magnetic phase 
is extended over the whole sample volume independently 
on the P value for all the investigated materials. The 
T-dependence of for the three samples at different 
values of P are reported in the main panels of Figs. [3j [4] 
and |Hj (LaCoPO, PrCoPO and LaCoAsO, respectively). 
Here, the To values are estimated as free parameters from 
the fitting to experimental ZF data according toEq. ([3]) 
where, in view of the results of ambient pressure mea- 
surements, (3 = 0.34 is kept as a fixed parameter for all 
the samples. 

By focussing on the results for LaCoPO first and by 
comparing the main panels of Figs. [2]and[3j it is imme- 
diate to realize that the qualitative effect of increasing P 
is extremely similar to what is induced by the full Pr/La 
substitution. In particular, a strong suppression of B^(0) 
together with a sizeable enhancement of Tc are obtained 
at the maximum applied P ~ 23.6 kbar. As shown by the 
dashed lines in the main panel of Fig. [3j the fitting func- 
tion reported in Eq. ^ still well reproduces the observed 
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FIG. 4: (Color online) Main panel: B M (T) at different P 
for PrCoPO. Dashed lines are best fits to experimental data 
according to Eq. [3] with f3 = 0.34 as fixed parameter. Close 
and open symbols for ambient P data refer to measurements 
performed at Dolly and at GPD with unloaded cell, respec- 
tively. Inset: Tc vs. P. The continuous line is a best-fit to 
data according to a linear function. 



FIG. 5: (Color online) Main panel: B M (T) at different P 
for LaCoAsO. Dashed lines are best fits to experimental data 
according to Eq. [3] with (3 = 0.34 as fixed parameter. Close 
and open symbols for ambient P data refer to measurements 
performed at Dolly and at GPD with unloaded cell, respec- 
tively. Inset: Tc vs. P. The continuous line is a best-fit to 
data according to a linear function. 



experimental data with f3 = 0.34 at all the P values. By 
more carefully investigating the experimental results re- 
ported in Fig. [3] one realizes that the suppression of is 
actually non-monotonic with increasing P (see also Fig. 
[8] later on). Remarkably, this sharp jump in the inter- 
nal magnetic field is not reflected in the T-dependence 
of Tq that steadily increases in a linear fashion across 
the whole experimental range (see the inset of Fig. [3]). 
Quantitative data relative to the linear increase of Tc as 
a function of P are reported in Tab. [TTJ 

As already shown in Fig. [2] for the measurements at 
ambient P, the T-dependence of B^ in PrCoPO is quali- 
tatively identical to what it is reported for LaCoPO. Re- 
markably, measurements performed at GPD (reported in 
Fig. [4]) confirm this qualitative result also for P values 
up to 23.4 kbar. As it is found in the case of LaCoPO, 
a linear increase of Tq as a function of P can be inferred 
(see the inset of Fig. [4] and Tab. [Tl]). However, differently 
from LaCoPO, the saturation value of the internal field 



TABLE II: Summarizing quantities of interest for the three 
investigated compounds after ZF-/i + SR measurements under 
P. Tc(0) represents the critical temperature at ambient P 
while ( dT c/dp) is the slope of the linear trends of Tc vs. P 
presented in the insets of Figs. |3j [4] and [5] 



Compound 


Tc(0) (K) 


(Vtc(0)) • (""fa/dP) (kbar" 1 ) 


LaCoAsO 


53.5 ± 1.0 


(2.35 ±0.15) x 10" 3 


LaCoPO 


33.2 ± 1.0 


(8.0 ±0.5) x 10~ 3 


PrCoPO 


48.0 ± 1.0 


(8.5 ±0.5) x 10" 3 



T>(0) is steadily enhanced upon increasing P, reaching 
5^(0) ~ 189 G at the maximum P value (see also Fig. [8] 
later on). 

A totally different phenomenology is detected for La- 
CoAsO up to the value P = 22.8 kbar. As it is clear 
from Fig. [5] the value -B M (0) — 2640 G does not display 
any dependence upon P within the experimental error 
while the linear increase of Tc is much less marked than 
in the case of LaCoPO and PrCoPO (see Tab. [TTf. It 
should be remarked that, as it is discussed in detail later 
in Sect. 



IV 



the sizeable difference of -B M (0) between the 
case of P- and As-based samples can be explained by the 
preferential occupation of different crystallographic sites 
by the muons in the two cases (close to LaO and Co(Pn) 
tri-layers, respectively). However, the function reported 
in Eq. ^ still yields to good fitting results for T>(T) 
data with j3 — 0.34 as a fixed parameter independently 
on the P value. Since both the materials with La as 
non-magnetic RE ion then share the same power-law-like 
trend, it is further confirmed that this must be entirely 
associated to the Co sublattice. It seems safe to deduce 
that this behaviour for B^ (T) is the fingerprint of static 
magnetism from Co in RECo(As,P)0 materials. Data at 
ambient P are in good agreement with what was previ- 
ously reported in the literature!^! 

In order to stress the strong qualitative analogies 
among the three investigated samples, the experimental 
data displaying B^ vs. T at all the different P values are 
plotted in Fig. [6] Here, a scaling procedure along both 
the T and £? M axes with, respectively, Tc and B^(0) val- 
ues proper of each data set allows us to clearly enlighten 
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FIG. 6: (Color online) vs. T for the three samples at all 
the investigated P values. Data clearly collapse on a single 
common power-law-like trend with j3 — 0.34 (continuous line) 
after a proper scaling of both T and _B M axes with Tfc and 
Bfj,(0) values, respectively. 



FIG. 8: (Color online) Main panel: B M vs. P at T = 5 K for 
both LaCoPO and PrCoPO. A biasing value Pb = 16.5 kbar 
was added to real P values for PrCoPO. The dashed line is 
a linear guide for the eye. Inset: values of fj,c deduced for 
LaCoPO and PrCoPO from the dc magnetization measure- 
ments presented in Fig. [7] The dashed line is a guide for the 
eye. 



a shared power-law-like trend with exponent j3 = 0.34 
common to all the samples and independent on P. 

As already comme nted above (and as further discussed 
in detail later in Sect 



IV) 



the relative variations of in 



ferromagnets cannot be ascribed only to variations in the 
value of the magnetic moment in the ordered phase and 
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FIG. 7: (Color online) Main panel: M vs. H at fixed T = 5 
K for LaCoPO at different P values and for PrCoPO at ambi- 
ent pressure after the subtraction of the linear paramagnetic 
term due to Pr 3+ ions. Inset: raw magnetization data for 
PrCoPO before the subtraction of the linear paramagnetic 
term evidenced by the dashed line. 



a further experimental study for such quantity should 
be independently performed. To this aim, the behaviour 
of /ico vs. P in LaCoPO was investigated in the low- 
P region by means of dc magnetometry. Measurements 
of dc magnetization (M) as a function of the magnetic 
field (H) were performed at T = 5 K up to P = 10 
kbar. The results are reported in the main panel of Fig. 
[7] The unique intrinsic contribution to M in LaCoPO 
is expected to come from the weakly-itinerant FM state 
associated with the Co sublattice and an ordered value 
of /ico C ° P ° = 0.3 ± 0.02 fi B can be in deed estimated, in 
good agreement with previous reports!- 17 ! 28 ! Remarkably, 
/1qq CoPO turns out to be independent on P within the ex- 
perimental error in the investigated P range. The same 
scenario holds for PrCoPO, as it is deduced from M vs. 
H measurements at ambient P (raw data are presented 
in the inset of Fig. [7]) . After the subtraction of the linear 
term enlightened in the inset of Fig. [7] by means of a 
dashed line and accounting for the paramagnetic contri- 
bution of Pr 3+ magnetic moments, one obtains the intrin- 
sic contribution of the Co sublattice reported in the main 
panel of Fig. [7] Results only show a slight reduction of 
the magnetic moment of Co in PrCoPO with respect to 
what was found for LaCoPO (fJ%? oPO = 0.275±0.025 // B 
should be compared with Aico C ° P ° = °- 3 ± °- 02 Mb)- 
These results clearly indicate that the dramatic drop in 
(P) described in the main panel of Fig. [3] cannot be 
ascribed to a sudden suppression of fico- 

A summary of the main experimental results for B^ vs. 
P in LaCoPO and PrCoPO at T = 5 K is reported in 
the main panel of Fig. [8j Data already reported in Fig. 
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[3] for LaCoPO have been complemented by two single 
measurements at T = 5 K and at low P values (~ 3 and 
~ 8 kbar). This allows us to confirm the reproducibil- 
ity of the dramatic reduction of and to characterize 
it as a sudden phenomenon occurring in a window with 
maximum width of few kbar around P ~ 5 kbar. Data 
for PrCoPO, on the other hand, have been plotted in the 
main panel of Fig. [8] after adding a phenomenological 
biasing pressure Pb = +16.5 kbar. It is interesting to 
notice that, after such shift along the P axis, a system- 
atic common enhancement of S„ with increasing P can 
be clearly discerned sharing the same slope (1.95 ± 0.05 
G/kbar) for the two different materials. Naively, one can 
think of Pb as an effective pressure equivalent to the in- 
ternal chemical pressure associated with the full Pr/La 
substitution. The inset of Fig. [8] summarizes the be- 
haviour of //Co vs. P as determined by dc magnetization 
measurements. It is clear how, in the range of P values 
up to 10 kbar, the value of fico is indeed constant and 
approximately equal to 0.3 /iB- 

Modifications of contributions to the local field other 
than the dipolar coupling among /i + and Co layers should 
then be responsible for the observed behaviour. This 
topic will be discussed in more detail in Sect. [TV] from 
a computational point of view. In particular, it will 
be shown how a sizeable effect of P in modifying the 
band structure and, accordingly, the transferred hyper- 
fine term to the \x + must be taken into consideration in 
order to explain the observed behaviour. 



IV. INSIGHTS INTO THE INTERNAL FIELD 
AT THE MUON SITE FROM FIRST-PRINCIPLE 
CALCULATIONS 

In this Section we will discuss the sharp suppression 
of Bfj,(P) occurring at P = 5 ± 2 kbar in LaCoPO (as 
presented in the main panel of Fig. [8]) by taking into 
consideration different possibilities explaining its origin. 
To this aim, one should first write down the internal field 
at the fi + site (whose position is identified by the 
vector Tp) in a non-magnetized FM material as a sum of 
several contributions. In particular, one has 

B, t = |B dip (r A1 )+B c (r Al ) + B L | (4) 

where Bdi p (r M ) is the dipolar field arising from the atoms 
within a sphere with diameter smaller than i? m( j {R m d be- 
ing the approximate linear dimension of a magnetic do- 
main), B c (r M ) is the contact hyperfine field and Bl the 
Lorentz field. By considering the value /Ltc — 0.3 /^b de- 
rived in LaCoPO (the superscript was dropped to the aim 
of clarity) by means of dc magnetization measurements, 
it is straightforwardly derived that |Bl| = ( 4 V 3 ) Msat — 
170 G. 

As previously discussed in Sect. |III[ no modification of 
HCo could be discerned in LaCoPO upon increasing the 
value of the external P. The drop of B^ measured at 



P = 5 ± 2 kbar may then only be related to the quanti- 
ties whose values could still be a function of P. namely 
r M and B c . In the former case, results would then im- 
ply that a change of the n + site is at the origin of the 
observed discontinuity. Such /i + jump could be possibly 
due either to the lattice undergoing a structural transi- 
tion or, more generally, to a smoother modification of 
the electrostatic energy landscape triggered by P induc- 
ing a change in the fi + preferential site. In this respect, 
it should be remarked that estimates of the structural 
parameters as a function of P by means of X-ray diffrac- 
tion measurements would be of the utmost importance 
in order to get further information. In the latter case, 
on the other hand, the sudden modification of the B c 
term could follow from the evolution of the structure of 
electronic bands upon increasing P. One could also con- 
sider other possibilities like, for instance, a P-induced 
change in the magnetic GS from e. g. FM to AFM, or 
a reorientation of the magnetic configuration in the Co 
sublattice upon gradually increasing P. Finally, and less 
interestingly, another possibility is that the experimental 
findings are the result of the local distortion introduced 
by /i + once thermalized inside the material. 

In order to get more insights into the origin of the 
sharp drop in B^ vs. P as revealed by /i + SR experiments 
on LaCoPO, ab-initio calculations were then carried out 
in order to check for the reliability of all the scenarios 
proposed above for the material under investigation. The 
main objects were the computation as a function of P of 
the crystallographic sites where fi + sit and of the energy- 
volume curves allowing one to derive the structural and 
magnetic GSs of the system. 



A. Computational details 

We used both the plane wave (PW) and the full po- 
tential linearized augmented plane waves (FP-LAPW) 
methods as implemented in the VASPSl and EffPSl pack- 
ages. The Perdew-Bruke-Ernzerhof (PBE) functional 
was used in order to evaluate the exchange-correlation 
potential.^ As for the PW approach, the electron-ion 
interaction was described by the projector augmented 
wave (PAW) pseudopotentials method.^ Electronic con- 
vergence was set up at 10~ 6 eV and sampling of the Bril- 
louin zone performed by the Monkhorst-Pack scheme^ 
on a 8 x 8 x 4 grid. A plane-wave cutoff of 600 eV and 
a Gaussian smearing of 0.02 eV was used throughout. 
FP-LAPW calculations were carried out using a basis 
set with RM? x max(||fc||) = 7.5, R™? being the sma u_ 
est muffin-tin (MT) radius inside the MT spheres, and 
( max = 8 for the angular momentum expansion in the 
MTs (for both the wave functions and the potential). 
The reciprocal space was sampled with the same grid 
used in the PW approach. A close agreement is observed 
between the results obtained by means of both the com- 
putational methods. 

Energy-volume curves were obtained from the PW- 



TABLE III: Crystallographic positions ol the interstitial sites A, B and C at two different values of P. The absolute values of 
the local magnetic field at A, B and C arising only from the Bdi P contribution in Eq. Q are reported for two different FM 
configurations of fico- The values of the minima of potential energy for the three sites are denoted by Vq while E% P represents 
the lower eigenvalue for the ZP motion of ji + into the i potential dip (see text). The values for the energy barriers between 
different sites are denoted as A iti i. All the energy values (Vq, E^p and A i 4 /) are conventionally referred to V A . 



Site 


Wyckoff position 


x 00, v ( a ), z ( c ) 


(G) 


BZ° ±C (G) 


VS (eV) 


£zp (eV) 


Ambient P (A AB = A B c = 1-42 eV, A A c = 1-06 eV) 


A 
B 
C 


2c 
2c 
4/ 


1/4 1/4 0.875 
1/4 1/4 0.42 

3/4 1/4 0.30 


330 
1100 


170 

550 



0.27 
0.56 


0.45 
0.63 
1.13 


P = 30 kbar (A AB = A BC = 1.36 eV, A AC = 1-08 eV) 


A 
B 

C 


2c 
2c 
4/ 


1/4 1/4 0.875 
1/4 1/4 0.42 

3/4 1/4 0.30 


330 
1100 


170 

550 



0.37 
0.62 


0.45 
0.73 
1.22 



based calculations by constant volume energy minimiza- 
tion. The convergence criteria for forces minimization 
was set to 5 x 10~ 3 eV/A and 10 -5 eV was used as 
threshold for self-consistent electronic cycles. The op- 
timized unit cell volume at ambient P for LaCoPO is 
133.21 A 3 with lattice parameters a = 3.966 A and 
c = 8.468 A. In accordance with previous findings,^ 
the density-functional theory (DFT) calculations repro- 
duce the experimental structural parameters with errors 
~ 1%. No anomalies in the energy- volume curves are ob- 
served within the explored P range. At the same time, 
the FM-ordered configuration is the GS for P values up 
to 100 kbar in agreement with experimental results.^ 
The results then suggest that both the crystal struc- 
ture and the FM configuration as GS are stable against 
the increase of P in LaCoPO at least in the investi- 
gated P range. The calculated value for the magnetic 
moment /^co = 0.57 fiB is in agreement with previous 
calculations. 17 This value is slightly higher than our ex- 
perimental estimate (ico — 0.3 and it is substantially 
unchanged for P < 40 kbar while, for higher pressures, 
the magnetic moment on Co atoms linearly decreases 
reaching 0.55 hb for P = 100 kbar. 



B. The interstitial fi sites and their stability 
against zero-point motion and external pressure 

The crystallographic sites where /i + stop after ther- 
malization processes may be identified to a first extent 
by calculating the GS electron density of the bulk mate- 
rial under investigation. The minima of the electrostatic 
potential consequently obtained are indeed found to pro- 
vide a correct estimate to this aimP^ Three inequivalent 
positions are computed for LaCoPO, whose absolute val- 
ues of the potential energy are denoted as Vq (i = A, B, 
C) where we conventionally set V A = eV. Site A (B) 
is located within the LaO (CoP) tri-layers while site C is 
aligned with O and Co in between the different tri-layers, 
as shown in Fig. [9b. Remarkably, we note that the elec- 



trostatic interaction favours site A unlike what is found 
in RECoAsO and REFeAsO where the inte rstitia l site 
close to the transition metal plane is favoredP^^ More 
detailed information about the crystallographic positions 
of the three sites and their evolution upon increasing P 
is reported in Tab. Ill and in Fig. [9] The local mag- 



III (assuming an 



netic fields at the /i + sites Bdi P arising from the dipolar 
contribution of /ico — 0.3 /ib magnetic moments were 
also computed and reported in Tab. 
undistorted lattice). 

Once thermalizcd, in view of the quite high zero-point 
(ZP) motion energy for the potential minima are not 
necessarily stable trapping sites. Whether the ZP energy 
is higher than the potential barrier surrounding a mini- 
mum, fi + can escape to lower-energy ones. To check if 
the minima yield stable trapping, we computed the ZP 
motion by solving the Schrodingcr equation for the fi + 
in the electrostatic potential. Then, we inspected the 
position of the probability maxima for the relative GS 
wavefunctions and characterized each fj, + site by the rel- 
ative lowest eigenvalue Ezp (see Tab. III). It should be 



remarked that a stable trapping site must have a node- 
less and well-localized GS. Calculations show that this is 
clearly the case for sites A and B, while for site C we 
found a delocalized wavefunction with a large amplitude 
around C (resonant state). Moreover for site C, at vari- 
ance with the findings for sites A and B, the lowest eigen- 
value E% P is higher than the energy barrier separating A 
and C. Under these conditions, all /i + stopping in C reach 
the energetically more favorable minimum A. We cannot 
exclude the presence of a self trapping mechanism which 
could, in principle, prevent fi + from reaching the lowest 
energy site A from minimum C. Anyway, the C site has 
a so high energy that the fraction of fi + stopping there 
should be negligible. From the above considerations, only 
A and B minima are eligible as stable /i + sites and C will 
be disregarded in the rest of the discussion. 

As already mentioned, DFT calculations does not evi- 
dence structural or magnetic transitions triggered by P. 
The energy minima and ZP energies of the interstitial 
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FIG. 9: (Color online) Isosurfaces of the electrostatic potential of LaCoPO for ambient P (upper) and P = 30 kbar (lower). For 
each pressure, the energies corresponding to the first eigenvalue for a muon in site A, B, and C are shown, i.e., V(r) iso = E% P 
for a) and d); V(r) iso = £§ P for b) and e); V(r) iso = E% P for c) and f). 



sites are slightly modified with increasing P up to 30 
kbar but variations are smooth and of the order of the 
estimated accuracy of the calculation. The same trend 
shown in Tab. Ill for the energy minima and the zero 



point motion is found also for higher pressures (P up 
to 100 kbar). Minima B and C, in particular, increase 
their energy with respect to site A while the ZP ener- 
gies of both the sites (with reference to the relative en- 
ergy minimum) do not change significantly. This allows 
to state that the occupation probability of those sites 
should remain practically unchanged upon increasing P 
or, at least, that site A turns out to be even more and 
more favourable from an energetic point of view in com- 



parison with B and C. In conclusion, no evidence for 
sudden changes in the population of interstitial sites can 
be derived from DFT calculations. 

After considering both the experimental and the com- 
putational results, the possibility of a sudden change of 
site occupation for fi + upon increasing P seems to be 
very unlikely, as described in Sect. |IV| From the exper- 
imental side, in fact, no sign of occupancy of different 
sites is present for any of the investigated samples (see 
the inset of Fig. [2]for representative raw data). An unre- 
alistically complete redistribution of the site occupation 
between the two inequivalent sites A and B would then 
be needed in order to explain the experimental results. 
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As already stressed above, moreover, no sign of struc- 
tural transitions upon increasing P for the investigated 
samples was evidenced by means of DFT computations. 
Due to the smoothness of the energy landscape coming 
out from calculations, then, it seems to be highly unlikely 
to attribute the sharp and dramatic drop in the internal 
field occurring at P ~ 5 kbar (see the main panel of Fig. 
[8l) to a modification in r„. 



C. Computation of the perturbation effect induced 
by the muon 

In this paragraph we will consider the local perturba- 
tion induced by the implanted muons. We studied the ef- 
fect of fi + induced lattice distortion on the electrostatic 
potential landscape. Normally, the dielectric screening 
in metals is so efficient that /i + will not cause significant 
lattice distortion. However, here we are dealing with a 
material that is a poor metal, and so it is important to 
get an estimate of any lattice distortion effect. Since we 
need to get just an estimate we can ignore the effect of ZP 
motion and treat the /z + as if it were a hydrogen inter- 
stitial impurity. We model the isolated impurity within 
the supercell approach by building a 64 atoms supercell 
up from our bulk structure. The supercell volume was 
kept constant while atomic positions were allowed to re- 
lax. /i + are positively charged particles but, at the same 
time, our supercell must have a neutral charge because 
the electronic screening in metals does neutralize the /z + 
charge. For the hydrogen impurity we chose those sites 
identified as absolute minima by the electrostatic poten- 
tial landscape. We then let the system evolve towards 
the GS allowing both electron rearrangement and lattice 
distortion. The final optimized position for the impurity 
represents the refined site for the and we found that 
the refined /i + position is in agreement with the one ob- 
tained by analyzing the electrostatic potential minima. 
Fully relaxed structures at ambient P, at P = 15 kbar 
and at P = 30 kbar show that both the /i + position in- 
side the cell and the distance between /i + and P ions 
varies less than 0.04 A. The phosphorous ion close to /i + 
is pushed closer to the Co plane by ~ 0.06 A and its 
four neighbouring Co atoms increase their magnetic mo- 
ment to 0.6 hb- Nevertheless, once more, no appreciable 
modification of the crystal structure and of the magnetic 
properties of the whole system (crystal and which 
could in principle justify the drop of the internal field 
observed around P ~ 5 kbar could be computed as a 
function of P. 
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FIG. 10: (Color online) Upper panels: energy bands of La- 
CoPO at ambient P and at 100 kbar (E F = 0). Black(red) 
color refer to majority (minority) spin bands. Lower panel: 
difference between the energy at Y and the Fermi energy Ep 
for the band crossing the Fermi energy. Inset: energy band 
dispersion for selected P values. 



kbar an unoccupied electron band shifts across the Fermi 
energy as shown in Fig. [lO] creating a large cylindri- 
cal Fermi surface owing to the flat dispersionless band 
along T-Z shown in figure. This modification in the band 
structure strongly suggests that the variation of the lo- 
cal field upon increasing P can be due to a change in the 
transferred hyperfine contact field B c (r„) reported in Eq. 
Q, strongly influenced by the conduction electrons It 
should be remarked that, according to what is reported 
in Tab. |III| for the in-plane orientation of the spins one 
has Bdip ( r Ai(A)) = — Bx and the muon in site A is thus 
subject only to the contribution of B c (r M ). 



V. DISCUSSION 



D. Pressure-induced effects on the energy bands 

In order to identify possible changes of the spin den- 
sity at the ix + site we computed the band structures to 
evaluate the contribution of the conduction electrons to 
the hyperfine field as a function of pressure. At P ~ 38 



Both chemical and external pressure have been shown 
to affect the ferromagnetic properties of the Co-based 
1111 oxy-pnictides. ZF-/i + SR experiments display that 
the magnetic ordering temperature T c increases as a func- 
tion of pressure. The internal field at the muon site is 
found to abruptly decrease around 5 kbar in LaCoPO. A 
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change of the same amount is found when lanthanum is 
substituted by smaller Pr ions. This demonstrates that 
chemical pressure mimics the hydrostatic pressure as can 
be deduced from Fig. [8] Claims of close analogies among 
the effect of chemical and external pressures were already 
reported concerning both 1111 and 122 compoundspSHS 
In order to understand the origin of the drastic reduc- 
tion observed for we have investigated the behav- 
ior of those quantities, implicitly included in Eq. Q, 
that mainly contribute to the precession frequency of the 
muon around the local field, namely the Co ordered mag- 
netic moment, fJ-co, and both the distance and the density 
of the spin at the /i + site as a function of P. 

Magnetization measurements have shown no significa- 
tive change of fico neither by changing chemical or exter- 
nal pressure (see the inset of Fig. f7|. DFT calculations 
have been performed in order to determine the muon sites 
and displays that only two lattice positions are stable, 
namely close to the CoP (site A) and to the REO (site 
B) layers (Fig. |9j), the latter being energetically favored. 
Since ZF-/i + SR display only one oscillating component, 
only one site is really occupied. The calculation of elec- 
trostatic potentials does not display any strong variation 
as a function of pressure, indicating that the /i + site can- 
not change. Since also fico is constant, we can conclude 
that neither the dipolar interaction or the Lorenz field 
of Eq. Q can be responsible for the suppression of the 
local field at the fi + site as a function of pressure. On the 
other hand, DFT band structure calculations show that 
a slight variation in the energy of a minority band occurs 
as a function of pressure, displaying that a crossing of 
the Fermi level takes place around 38 kbar, as shown in 
Fig. [lOj This computational outcome suggests a sudden 
change of the hyperfme contribution in Eq. Q due to 
the occupation of this minority spin band. A notable 
discrepancy between the critical P value for the drop 
of the internal field from experimental (~ 5 kbar) and 
DFT (~ 38 kbar) results is found. This difference may 
be understood by considering that, first, the band struc- 
ture is extremely sensitive to the position of the P ions. 
The shift of the position of P ions towards the Co plane 
introduced by the incoming /i + discussed above is not 
taken into account when inspecting the band structure 
and it may indeed favour the population of the unoccu- 
pied valence band at lower P values. In addition, in this 
delicate scenario, even small non-hydrostatic effects not 
taken into account in the calculations can shift the criti- 
cal pressure. Moreover, the DFT approach for the first- 
principles analysis may underestimate P effects as a con- 
sequence of the overestimated magnetic moment which 
in turn is related to the height of phosphorous ionsP^l 
Unfortunately, as already recalled above, no experimen- 
tal structural refinements as a function of P are available 
at the moment. Thus, we cannot quantify the mismatch 
between DFT obtained and real crystal structures as a 
function of pressure. 

Other hypotheses can be made in order to explain 
the observed phenomenology by referring to the intrin- 



sic magnetic properties of the materials, whose changes 
upon increasing P can be directly reflected in relative 
modifications into one (or more) terms appearing in Eq. 
@. One can then take into consideration the possibility 
of a P-triggered reorienting of the FM configuration. The 
values for the modification of the overall local magnetic 
field at the muon site according to this hypothesis are 
provided in Tab. |III| Besides a possibly satisfying expla- 
nation of experimental data from a quantitative point of 
view, this hypothesis should be considered as unrealistic 
mainly due to the low absolute value of spin-orbit inter- 
action for the Co ions. The different energies of the GSs 
were analyzed by performing non-collincar spin-polarized 
calculations for different P and easy magnetization axes 
orientations. The results for different spin configurations 
are found to be nearly degenerate so that, in conclusion, 
the mean-field approach does not provide a clue to sup- 
port this alternative. 



VI. CONCLUSIONS 

We performed detailed measurements by means of 
muon spin spectroscopy on LaCoPO, PrCoPO and La- 
CoAsO under applied hydrostatic pressure. The localized 
electrons on the ionic shells of Pr 3+ ions do not affect at 
all the local static features of magnetism as detected by 
muons in spite of the high value of the magnetic mo- 
ment expected for the free Pr 3+ ion. Phosphorous-based 
compounds turn out to be much more sensitive than 
LaCoAsO to the application of pressure. In particular, 
the critical transition temperature to the ferromagnetic 
phase in LaCoPO is sizeably increased both by chemical 
and external pressures (the former being triggered by the 
full Pr 3+ /La 3+ substitution), while for LaCoAsO the ef- 
fect is also present but by far less marked. The increase 
of both kinds of pressure, moreover, dramatically sup- 
press the local magnetic field at the muon site leaving the 
magnetic moment per Co ion substantially unchanged. 
Density-functional theory calculations of the band struc- 
ture suggest that this the change of the muon field is 
due to a subtle change of the hyperfme coupling driven 
by a minority band crossing the Fermi surface. Results 
clearly evidence how crucial the role of pressure is on the 
properties of the ferromagnetic phase and, at the same 
time, how similar the effects of chemical and external 
hydrostatic pressure is in 1111 oxy-pnictide compounds. 
A computational description of LaCoPO performed by 
means of ab-initio DFT calculations supports these find- 
ings and suggests that chemical and external pressures 
both trigger a qualitative change in the electronic band 
structure of the investigated compounds. 



VII. ACKNOWLEDGEMENTS 

G. Prando acknowledges support from the Leibniz- 
Deutscher Akademischer Austauschdienst (DAAD) Post- 



12 



Doc Fellowship Program. G. Profeta acknowledges sup- 
port from the FP7 European project SUPER-IRON 
(grant agreement No. 283204), by a CINECA-HPC IS- 
CRA grant and by an HPC grant at CASPUR. F. B. ac- 
knowledges support from CASPUR under the Standard 
HPC Grant 2012 and from the FP7 European project 
SUPER-IRON (grant agreement No. 283204). E. M. 



B. and H.-J. G. acknowledge support from Deutsche 
Forschungsgemeinschaft (DFG) through SPP1458 (grant 
No. GR3330/2). P. C. and S. S. acknowledge support 
from Fondazione Cariplo (research grant No. 2011-0266). 
The structural models in the figures were rendered using 
the VESTA packageP 



* E-mail: g.prando@ifw-dresden.de 

1 Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. 
Am. Chem. Soc. 130, 3296 (2008) 

2 D. C. Johnston, Adv. Phys. 59, 803 (2010) 

3 Z.-A. Ren, W. Lu, J. Yang, W. Yi, X.-L. Shen, Z.-C. Li, G- 

C. Che, X.-L. Dong, L.-L. Sun, F. Zhou, and Z.-X. Zhao, 
Chin. Phys. Lett. 25, 2215 (2008) 

4 G. Prando, A. Lascialfari, A. Rigamonti, L. Romano, S. 
Sanna, M. Putti, and M. Tropeano, Phys. Rev. B 84, 
064507 (2011) 

5 A. S. Sefat, A. Huq, M. A. McGuire, R. Jin, B. C. Sales, 

D. Mandrus, L. M. D. Cranswick, P. W. Stephens, and K. 
H. Stone, Phys. Rev. B 78, 104505 (2008) 

6 A. Marcinkova, D. A. M. Grist, I. Margiolaki, T. C. 
Hansen, S. Margadonna, and J.-W. G. Bos, Phys. Rev. 
B 81, 064511 (2010) 

7 V. P. S. Awana, A. Pal, A. Vajpayee, R. S. Meena, H. 
Kishan, M. Husain, R. Zeng, S. Yu, K. Yamaura, and E. 
Takayama-Muromachi, J. Appl. Phys. 107, 09E146 (2010) 

8 T. Shang, L. Yang, Y. Chen, J. L. Zhang, L. Jiao, J. Chen, 
J. Dai, and H. Q. Yuan, |arXiv:1212.022i"1 (2012) 

9 T. Shang, L. Yang, N. Cornell, F. Ronning, Y. Chen, L. 
Jiao, J. Chen, A. Howard, J. D. Thompson, A. Zakhidov, 
M. B. Salamon, and H. Q. Yuan, |arXiv:1212.0222| (2012) 

10 A. S. Sefat, R. Jin, M. A. McGuire, B. C. Sales, D. J. Singh, 
and D. Mandrus, Phys. Rev. Lett. 101, 117004 (2008) 

11 A. Yamamoto, J. Jaroszynski, C. Tarantini, L. Balicas, 
J. Jiang, A. Gurevich, D. C. Larbalestier, R. Jin, A. S. 
Sefat, M. A. McGuire, B. C. Sales, D. K. Christen, and D. 
Mandrus, Appl. Phys. Lett. 94, 062511 (2009) 

12 G. Prando, R. Giraud, M. Abdel-Hafiez, S. Aswartham, 
A. U. B . Wolter, S. Wurmehl, and B. Biichner, 
|arXiv:1207.2457| (2012) 

13 H. Maeter, J. E. Hamann-Borrero, T. Goltz, J. Spehling, 
A. Kwadrin, A. Kondrat, L. Veyrat, G. Lang, H.-J. Grafe, 
C. Hess, G. Behr, B. Biichner, H. Luetkens, C. Baines, A. 
Amato, N. Leps, R. Klingeler, R. Feyerherm, D. Argyriou, 
and H.-H. Klauss, |arXiv: 1210.6959] (2012) 

14 G. Prando, P. Carretta, A. Rigamonti, S. Sanna, A. Palen- 
zona, M. Putti, and M. Tropeano, Phys. Rev. B 81, 
100508(R) (2010) 

15 H. Maeter, H. Luetkens, Y. G. Pashkevich, A. Kwadrin, R. 
Khasanov, A. Amato, A. A. Gusev, K. V. Lamonova, D. A. 
Chervinskii, R. Klingeler, C. Hess, G. Behr, B. Biichner, 
and H.-H. Klauss, Phys. Rev. B 80, 094524 (2009) 

16 R. De Renzi, P. Bonfa, M. Mazzani, S. Sanna, G. Prando, 
P. Carretta, R. Khasanov, A. Amato, H. Luetkens, M. Ben- 
dele, F. Bernardini, S. Massidda, A. Palenzona, M. Tro- 
peano, and M. Vignolo, Supercond. Sci. Technol. B 25, 
084009 (2012) 

17 H. Yanagi, R. Kawamura, T. Kamiya, Y. Kamihara, M. 
Hirano, T. Nakamura, H. Osawa, and H. Hosono, Phys. 



Rev. B 77, 224431 (2008) 

18 H. Ohta and K. Yoshimura, Phys. Rev. B 79, 184407 
(2009) 

19 H. Ohta and K. Yoshimura, Phys. Rev. B 80, 184409 
(2009) 

20 H. Ohta, C. Michioka, A. Matsuo, K. Kindo, and K. 
Yoshimura, Phys. Rev. B 82, 054421 (2010) 

21 H. Ohta, C. Michioka, and K. Yoshimura, Journ. Phys. 
Soc. Jpn. 79, 054703 (2010) 

22 R. Sarkar, A. Jesche, C. Krellner, M. Baenitz, C. Geibel, 
C. Mazumdar, and A. Poddar, Phys. Rev. B 82, 054423 
(2010) 

23 V. P. S. Awana, I. Nowik, A. Pal, K. Yamaura, E. 
Takayama-Muromachi, and I. Felner, Phys. Rev. B 81, 
212501 (2010) 

24 A. Pal, M. Tropeano, S. D. Kaushik, M. Hussain, H. Kis- 
han, and V. P. S. Awana, J. Appl. Phys. 109, 07E121 
(2011) 

25 H. Ohta, C. Michioka, and K. Yoshimura, Phys. Rev. B 
84, 134411 (2011) 

26 J. Sugiyama, M. Mansson, O. Ofer, K. Kamazawa, M. 
Harada, D. Andreica, A. Amato, J. H. Brewer, E. J. 
Ansaldo, H. Ohta, C. Michioka, and K. Yoshimura, Phys. 
Rev. B 84, 184421 (2011) 

27 C. Krellner, U. Burkhardt, and C. Geibel, Physica B 404, 
3206 (2009) 

28 A. Pal, S. S. Mehdi, M. Husain, B. Gahtori, and V. P. S. 
Awana, J. Appl. Phys. 110, 103913 (2011) 

29 M. Majumder, K. Ghoshray, A. Ghoshray, B. Bandyopad- 
hyay, B. Pahari, and S. Banerjee, Phys. Rev. B 80, 212402 
(2009) 

30 H. Sugawara, K. Ishida, Y. Nakai, H. Yanagi, T. Kamiya, 
Y. Kamihara, M. Hirano, and H. Hosono, Journ. Phys. 
Soc. Jpn. 78, 113705 (2009) 

31 M. Majumder, K. Ghoshray, A. Ghoshray, B. Bandyopad- 
hyay, and M. Ghosh, Phys. Rev. B 82, 054422 (2010) 

32 M. Majumder, K. Ghoshray, A. Ghoshray, A. Pal, and V. 
P. S. Awana, J. Phys. Soc. Jpn. 81, 054702 (2012) 

33 E. M. Briining, C. Krellner, M. Baenitz, A. Jesche, F. 
Steglich, and C. Geibel, Phys. Rev. Lett. 101, 117206 
(2008) 

34 Y. Luo, Y. Li, S. Jiang, J. Dai, G. Cao, and Z.-A. Xu, 
Phys. Rev. B 81, 134422 (2010) 

35 R. Sarkar, M. Baenitz, A. Jesche, C. Geibel and F. Steglich, 
Journ. Phys. Cond. Matt. 24, 135602 (2012) 

36 A. Jesche, T. Forster, J. Spehling, M. Nicklas, M. de Souza, 
R. Gumeniuk, H. Luetkens, T. Goltz, C. Krellner, M. Lang, 
J. Sichelschmidt, H.-H. Klauss, and C. Geibel, Phys. Rev. 
B 86, 020501(R) (2012) 

37 N. Colonna, G. Profeta, and A. Continenza, Phys. Rev. B 
83, 224526 (2011) 

38 S. A. J. Kimber, A. Kreyssig, Y.-Z. Zhang, H. O. Jeschke, 



13 



R. Valentf, F. Yokaichiya, E. Colombier, J. Yan, T. C. 
Hansen, T. Chatterji, R. J. McQueeney, P. C. Canfield, 
A. I. Goldman, and D. N. Argyriou, Nature Mater. 8, 471 
(2009) 

39 G. Prando, S. Sanna, G. Lamura, T. Shiroka, M. Tropeano, 
A. Palenzona, H.-J. Grafe, B. Biichner, P. Carretta, and 
R. De Renzi, |arXiv:1209.6512| (2012) 

40 E. Gati, S. Kohler, D. Guterding, B. Wolf, S. Knoner, 
S. Ran, S. L. Bud'ko, P. C. Canfield, and M. Lang, 
|arXiv:1210.5398| (2012) 

41 F. Nitsche, A. Jesche, E. Hieckmann, T. Doert, and M. 
Ruck, Phys. Rev. B 82, 134514 (2010) 

42 A. Yaouanc and P. Dalmas de Reotier, Muon Spin Ro- 
tation, Relaxation, and Resonance: Applications to Con- 
densed Matter, Oxford University Press, Oxford (2011) 

43 R. Khasanov, S. Sanna, G. Prando, Z. Shermadini, M. 
Bendele, A. Amato, P. Carretta, R. De Renzi, J. Karpinski, 
S. Katrych, H. Luetkens, and N. D. Zhigadlo, Phys. Rev. 
B 84, 100501(R) (2011) 

44 W. J. Duncan, O. P. Welzel, C. Harrison, X. F. Wang, X. 
H. Chen, F. M. Grosche, and P. G. Niklowitz, J. Phys.: 



Cond. Mat. 22, 052201 (2010) 

V. H. Tran, P. Rogl, P. Dalmas de Reotier and A. Yaouanc, 
Phys. Rev. B 83, 144417 (2011) 

S. Barth, E. Albert, G. Heiduk, A. Moslang, A. Weidinger, 
E. Recknagel and K. H. J. Buschow, Phys. Rev. B 33, 430 
(1986) 

G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 

(1996) ; 

Elk code (version 1.4.18), http://elk.sourceforge net 



J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 
77, 3865 (1996) 

P. E. Blochl, Phys. Rev. B 50, 17953 (1994) 

H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 

(1976) 

C. P. Slichter, Principles of Magnetic Resonance, Springer- 
Verlag Berlin (1990) 

Y.-Z. Zhang, H. C. Kandpal, I. Opahle, H. O. Jeschke, and 

R. Valentf, Phys. Rev. B 80, 094530 (2009) 

K. Momma, and F. Izumi, J. Appl. Cryst. 44, 1272 (2011) 



